Quantum inelastic conductance through molecular wires 
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We calculate non-perturbatively the inelastic effects on the conductance through a conjugated 
molecular wire-metal heterojunction, including realistic electron-phonon coupling. We show that at 
sub-band-gap energies the current is dominated by quantum coherent transport of virtual polarons 
through the molecule. In this regime, the tunneling current is strongly increased relative to the 
case of elastic scattering. It is essential to describe the full quantum coherence of the polaron 
formation and transport in order to obtain correct physics. Our results are generally applicable to 
one-dimensional atomic or molecular wires. 



PACS numbers: 85.30.Vw, 73.50.-h, 73.40. Gk, 73.61.Ph 



i 

-a 
c 

o 
o 



(N 
> 

o 

m 
m 
O 

&\ 
On 
-i— > 



i 

C 

o 
o 



X 



Truly one-dimensional conducting structures, where 
the confinement lengths for electrons in two directions 
are of the order of atomic diameters, are currently the 
subject of much experimental work. Examples include 
fullerene nanotubes H , organic molecules bonded to 
metallic electrodes g]7 and dangling-bond (DB) wires 
created by scanning tunneling microscope lithography 
of the H-saturated Si(001) surface ||[|]. The conduc- 
tance properties of such structures are of great impor- 
tance when one considers their potential as atomic or 
molecular scale electronic devices. There have been many 
recent theoretical investigations of electron transport in 
these systems: for example, Joachim et al. studied Xe 
atom 'wires' and through organic molecules || us- 
ing the 'elastic scattering quantum chemistry' approach; 
idealised one-dimensional atomic wires sandwiched be- 
tween jellium electrodes have been studied using den- 
sity functional theory Jt| and a recursion transfer matrix 
method ||, while Datta and collaborators studied or- 
ganic molecules on metal electrodes [|| . 

All the above theoretical work was done within the 
approximation of elastic transport, where the electron- 
phonon (e-ph) coupling plays no role. But there are good 
reasons to doubt the validity of this approximation in 
one-dimensional atomic or molecular scale systems. In 
small systems, the coupling between electron and other 
excitations is enhanced. Furthermore, a one-dimensional 
metal is generally unstable towards a Peierls distortion 
pfj| . Once such a distortion has occurred and produced 
a band gap, charges added into the system tend to self- 
localise and cause distortions of the system which lower 
the band gap. Such polaronic phenomena have been 
studied in conducting polymers for decades Jill ; elec- 
tron transport usually proceeds via tunneling into po- 
laron states arising from lattice fluctuations |12| . 

Motivated by this physics and the possibility of mea- 
suring the transport through one-dimensional atomic and 
molecular wires, we report in this letter the first calcu- 
lations of the inelastic electronic transport through an 
atomic-scale wire including realistic e-ph coupling. In 



our calculations, the quantum coherence of the states 
is fully retained, i.e., no adiabatic separation between 
the electronic and phonon degrees of freedom is made. 
We have chosen to study a conjugated molecular chain 
as an example since the e-ph coupling in such materials 
is relatively well understood. However the same physi- 
cal principles (and therefore the qualitative features of 
our results) are widely applicable to all one-dimensional 
atomic and molecular wires. 

For energies within the gap, the coherent transport 
is dominated by tunneling. Here we study the equiva- 
lent of the polaronic phenomena in conducting polymers 
for virtual (i.e. tunneling) electron injection. This has 
previously been included in very few transport calcula- 
tions: the elastic conductivity was found for a conjugated 
molecule containing static solitons [jl3| and the classical 
response of a conjugated oligomer to an injected elec- 
tron wavepacket has been considered [ fl4j |. However, the 
quantum coherence between the electron and the lattice 
was not retained. Such coherence is kept in our approach 
since it is essential to treat the tunneling of objects such 
as polarons or solitons. 

Our approach for studying the coupled e-ph system 
is inspired by the Su-Schrieffer-Heeger (SSH) Hamilto- 
nian fll|| , in which a tight-binding treatment for the 
7r-electrons is combined with a classical ball-and-spring 
model for a one-dimensional atomic (CH)^ chain. In our 
model, the phonons are treated on a quantum level, the e- 
ph coupling is linear, and the electronic part is expressed 
in terms of the one-electron eigenstates, labelled by n. 
The Hamiltonian is 



H = J2 £ «< 



c„c„ + 



q,n,m 



(i) 



where c\ (c„) creates (annihilates) an electron in the n-th 
electronic state with energy e„ and a J (a q ) creates (anni- 
hilates) a phonon in mode q. This is a general Hamilto- 
nian, describing coupled e-ph systems going beyond the 
Holstein model. 
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In this work, the electronic energies e„, the phonon 
frequencies uj q and the coupling matrix elements 7 gnm 
are obtained from the SSH Hamiltonian for an iso- 
lated molecular chain The one-electron eigenstates 
and energies e n are calculated self-consistently with the 
atomic distortions for a given number of atoms iV a and 7r- 
electrons N c . The ground state of the neutral (N c = N a ) 
dimerised chain is our reference system, from which the 
phonon modes and frequencies are obtained by pertur- 
bation theory The energy change, taken to second 
order in the atomic displacements Xi , gives the dynamical 
matrix from which phonon eigenmodes V q (i) and frequen- 
cies are obtained. We have checked that this second-order 
expansion of the elastic energy is sufficient to describe the 
formation of a static polaron in the chain. The e-ph cou- 
pling matrix elements are calculated from the variation 
of the electronic Hamiltonian due to phonon displace- 
ments. They are transformed from a real-space represen- 
tation into the electron eigenstate representation to get 
the 

Iqnm elements in Eq. (|l|). 
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FIG. 1. Dimerisation pattern di = (— l) l (xi+i — 2xi+Xi-i) 
(in A) induced by an additional electron in the isolated wire 
(TVa = 40, iV c = 41): classical SSH model (dashed line); 
half electronic spectrum and all classical phonons (solid line); 
N p h = 6 (lowest optical modes) classical phonons (circles), 
N p h = 6 (same modes) quantum phonons with n™^ x = 6 
(dot-dashed line). 

We choose the eigenstate representation because we 
wish to consider the action of H on the (_/V a + l)-electron 
(or (N a — l)-electron) Hilbert space obtained by adding 
(or removing) one electron to (or from) the ground state 
of the neutral chain. It is most straightforward to project 
out the addition of an electron (or hole) to the already 
occupied (or empty) states if we work entirely with the 
eigenstates. Then the n,m sums in Eq.(^) run over 
the occupied valence band states (for hole transport) or 
over the unoccupied conduction band states (for electron 
transport). We note that it is essential to keep the N a 
electrons of the reference system implicitly in the calcu- 
lation since they drive the Peierls distortion in the sys- 
tem. We justify the inclusion of only (A^ a ± l)-electron 
states in the transport by noting that the mean time be- 



tween electron (or hole) passages (« 10~ 7 s for a current 
of 1 pA) is much longer than a typical residence time in 
the molecule (~ 10 -15 s for a molecule strongly coupled 
to the electrodes [0). Furthermore, we have verified 
that the static lattice distortions due to charging are well 
reproduced by considering only one half of the electronic 
spectrum (Fig. |l|). 

The degree of dimerisation is conveniently represented 
by the staggered difference di between adjacent bond 
lengths. Figure |l| shows the di pattern induced when one 
additional electron is injected into the neutral isolated 
chain. As expected, we find a reduction of dimerisation, 
corresponding to a polaron-like defect in the middle of 
the chain. For short chains (A r a =40 for example), the 
difference of the dimerisation calculated with (a) the full 
electronic spectrum and an exact classical treatment of 
the lattice (SSH model), and (b) half of the electronic 
spectrum and a harmonic treatment of the lattice, does 
not exceed 7% . We also checked that only the longest 
wavelength optical phonon modes Jl^] contribute to the 
polaron distortion. 

The basis states for the Hamiltonian Eq.([|) are ex- 
pressed as \n, {n q }) where n indexes the one-electron 
eigenstates and {n q } represents the set of occupation 
numbers for the different phonon modes q. The elastic 
states \n, {0}) have no excited phonon modes and the in- 
elastic states \n, {n q }) have phonon modes q excited with 
n q = 1, 2, 3, .., n.™£ X . In all the calculations, the dynami- 
cal (i.e. energy dependent) correlation between the elec- 
tron and the phonons is fully kept. Even for a relatively 
short wire the size of the complete basis set quickly be- 
comes intractable for reasonable numerical calculations. 
However using a limited but sufficient number of phonon 
modes to describe the polaron formation permits us to 
treat a range of chain lengths with reasonable comput- 
ing times. Furthermore, our results show that only a 
few excitations in each phonon mode are needed to give 
converged results. 

Knowing e„, uj q and j qnm for a realistic description 
of the wire, we then calculate the transport properties 
through the electrode/wire/electrode junction. The con- 
jugated chain is connected to two semi-infinite ID metal- 
lic electrodes of band width 6 eV via the use of embedding 
potentials El and E^ [|o| . The propagation of an elec- 
tron or hole through the junction is treated as a scatter- 
ing problem with many channels [pl|-p3"| . Full scattering 
boundary conditions are applied for an incident electron 
or hole plane wave with the wire in a given vibrational 
state. For the results given in this paper, the chain config- 
uration before scattering is the ground state (no phonons 
present). The solution of the scattering problem can be 
obtained by solving the following linear system 

[E — H — Z L (E) - E R (E)] \a(E)) = \s(E)) , (2) 

where \s(E)) is the source term, representing the injected 
carrier at energy E. The ket \a(E)) is the wave function 
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of the scattering state; it is represented as an expansion 
over the basis set \n, {n q }). We can derive the transmis- 
sion and reflection coefficients for the different channels, 
the currents and the average atomic displacements from 
the solution (22|j2|] of Eq.(g) obtained by a conjugate 
gradient technique. 

To check the validity of our quantum approach, we 
first calculate the mean atomic displacements (xi) in- 
duced by adding an electron to an isolated chain in the 
ground state ]2^ j. The corresponding dimerisation pat- 
tern, also shown in Fig.[j], reproduces the results obtained 
with classical phonons. The dip in the dimerisation is 
slightly shallower in the fully quantum case; this may 
result from quantum delocalisation of the polaron state. 
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FIG. 2. Differential conductance o~(E) (in units of 2e 2 /h) 
through the wire (iV a = 40). cr(E) with inelastic scattering 
is represented by the symbols corresponding to different iVph 
optical modes and n™c"- The thin solid line is a guide for 
the eye. The elastic conductance through the static neutral 
chain (solid line) and charged chain containing a static po- 
laron (dashed line) is also shown (E — is at mid-band-gap). 

We now turn to calculations on a current-carrying wire 
between two electrodes. Figure || shows the differential 
conductance cr(E) of the junction. Calculations are per- 
formed with different numbers of phonon modes and with 
different maximum occupation numbers ri^?. We also 
calculate the elastic conductance through a static neutral 
chain and a charged chain containing a static polaron in 
the centre. One important result of our approach is that 
the inclusion of inelastic scattering strongly enhances the 
current in the tunneling regime compared to the case of 
elastic scattering. The transport becomes increasingly 
dominated by resonant tunneling through virtual polaron 
states as the injection energy approaches the top of the 
band gap. The first resonance in a(E) occurs at an in- 
jection energy very close to the charging energy i? c h of 
the isolated molecular chain (for N a — 40, we calculate 
Ech = 0.53 eV). As expected, for any injection energy 
within the tunneling gap, the current decreases exponen- 
tially as the length of the wire increases (Fig. |3|) for both 
elastic and inelastic calculations. The absolute value of 



(j(E) is much larger when the e-ph coupling is included, 
and the apparent band gap (related to the slope of the 
curve) is smaller. Such effects become more important 
for longer molecular wires (N a > 100). For shorter wires 
(A a < 10), both calculations seem to converge; the tun- 
neling process is too fast for the lattice to respond signif- 
icantly to the presence of the electron, so the effects of 
e-ph coupling are less pronounced H . 
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FIG. 3. Wire length dependence of u(E) (in units of 
2e 2 /h) for an injection energy E — 0.3 eV inside the tunnel- 
ing gap. The inelastic calculations are obtained for N p h = 4 
lowest optical modes, n™^ x = 2 (circles) and n™^ x = 3 (dia- 
monds). The elastic conductance through the static neutral 
chain is also shown (squares). 

Finally, we present in Figure the quantum averaged 
induced atomic displacement x^ol an atomic site j when 
the tunneling electron is found on site i j26|. These 
results are characteristic of phonon fluctuation-assisted 
tunneling via the dynamical formation of a virtual po- 
laron inside the chain. 
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FIG. 4. Induced atomic displacements (in A) due 
to a tunneling electron (injection energy at mid-gap) for 
a jV a = 100 chain. The curves are shifted vertically 
(by 0.02) for clarity. The electron-projection Pi sites are 
1=1,11,31,41,51,71,91 for the top to the bottom curve respec- 
tively. Calculations are performed for iVph = 8 optical modes 



The bond alternation, and hence the band gap, are 
locally suppressed around the tunneling electron. This 
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distortion accompanies the electron through the chain, 
leaving a 'wake' of distortion behind it. Such a persis- 
tent response of the lattice is not surprising in the ab- 
sence of damping p^|. The displacements are slightly 
less than those for an isolated chain (not shown) because 
the lattice cannot respond fully to the tunneling electron. 
This effect manifests itself in the calculated conductance 
(Fig. |2|); when tunneling, the current in the full quantum 
coherent calculations is less than the current that flows 
through a static chain where the polaron is 'frozen in'. 
This is because (a) the effective distortion in the quan- 
tum calculations is slightly smaller, and (b) the frozen-in 
polaron already exists so no elastic energy is needed to 
create it. Only fully quantum coherent calculations give 
qualitatively correct physics; it cannot be obtained from 
elastic scattering calculations through static chains (with 
or without disorder) or using a classical lattice. 

To conclude, we have presented the first calculation 
of inelastic effects on the conduction through molecular 
wires by explicitly considering realistic electron-phonon 
coupling. We have shown that the quantum coherent 
dynamical formation and transport of virtual polarons 
inside the molecular wire is crucial to describe the elec- 
tronic transmission properties of such a system, and 
leads to a strongly increased current in the tunneling 
regime. This is the small-bias regime in which molecular 
wires will normally operate. Therefore, whenever tunnel- 
ing conduction occurs in these or other one-dimensional 
atomic-scale wires [ p7[ and most probably in conducting 
polymer wires jl2| , it will be by the mechanism identified 
here. 
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